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?H ; Abstract 



We consider the exact effective superpotential of A/" = 1 U{Nc) super Yang-Mills theory 
^ I with Nf massive flavors an additional adjoint Higgs field. We use the proposal of Dijkgraaf 
and Vafa to calculate the superpotential in terms of a matrix model with a large number 
of flavors. We do this by gauging the flavor symmetry and forcing this sector in a classical 
vacuum. This gives rise to a 2-matrix model of ADE type A2, and large flavors. This 
approach allows us to add an arbitrary polynomial tree level superpotential for the Higgs 
field, and use strict large methods in the matrix model. 



1. Introduction 



Recently, Dijkgraaf and Vafa (DV) proposed an exact calculation of superpotentials in A/" = 1 
super Yang-Mills theories using matrix models P3I21E1- Later the derivations of the relation 
to the matrix models was shown directly from the field theory in j?l Ej. The DV proposal 
has initiated much research, extending the proposal to include flavors, other gauge groups, 
gravitational corrections, and relation to Seiberg Witten curves and field theory calculations 
of the effective superpotential 

The basic proposal of DV considers an U(A''c) gauge theory with an adjoint Higgs field 
$, which has a tree level superpotential Tr W{^). The solution is given in terms of a matrix 
model with of a M x M matrix $ with action j- Tr W{^). The effective superpotential for 
the glueball field S = {XX) can be calculated as the derivative of the planar contribution to 
the free energy in the matrix model 

WeffiS) = iV,^, (1) 

identifying S with the 't Hooft coupling QgM of the matrix model. As only the planar 
diagrams contribute, this can be calculated by large techniques. 

In the approach to the addition of flavors, there are two proposals. First in [Hj it is 
conjectures that one has to add the planar free energy with 1 quark-loop to the above 
effective action. This proposal has been used in most of the other approaches to the use of 
matrix models for super Yang-Mills with flavors, and gives agreement with previously known 
results. The proposal of jZj amounts to taking also the number of flavors Mf in the matrix 
model to infinity, but holding the ratio Mf/M fixed. One can easily convince oneselve that 
this gives contributions from all quark loops, so the proposals seem to be incompatible. One 
can argue that there should only contributions from planar diagrams with at most one quark 
loop to the effective superpotentials. This is a generalization of the arguments used in pi IH] 
to show that the effective superpotential has only contributions from the planar diagrams 
and is linear in A^^- What makes the proposal of jH] harder to work with is the fact that one 
needs 1/A^ corrections. These could in principle be calculated using techniques developed in 
[131113]. In most other papers on the subject only quadratic tree level superpotentials were 
considered, for which the planar diagrams can be summed explicitly. 

In this paper we will show that one can add flavors to DV staying purely in a large A^ 
limit. As in j7j we consider a matrix model with a large number of flavors Mf, but we will 
not constrain their ratio. This would add an extra large A^ parameter Sf = QsMf, which 
we will take to zero at the end of the calculation. The proposed formula for the effective 



1 



superpotential is then 

»K„(S)^A'.§ + A',||. (2) 

A simple counting of the dependence of the contributions on S and Sf shows that there 
are only contributions of the planar diagrams. Taking Sf = will then give only the 
contributions with at most one quark loop. Hence our proposal reproduces jU], but in a large 
context. 

Another way to study this system using only the original DV proposal can be found by 
gauging the flavor symmetry. For this we consider a SU(A'c) x SU(A^/) super Yang-Mills with 
bifundamentals and adjoint Higgs fields for each of the groups. The effective superpotential 
for this model can be calculated by an ADE type 2-matrix model as in [3 ES] . The form 
of the effective superpotential has roughly the form Q above. The SU(A''c) theory can be 
extracted by putting the SU(A^/) in a classical vacuum. Note that in this system Sf = QsMf 
becomes the gluino condensate in the SU(A^j) factor. The SU(A^/) Higgs field then will play 
the role the mass matrix for the quarks. This can be accomplished by adding and putting 
the SU(A^/) glueball field Sf to zero. As Sf = QsMf in the matrix model this corresponds 
nicely to what we had above. This system can be understood in terms of string theory by 
wrapping D5-branes on two different 2-spheres intersecting at a point, where one of the two 
spheres is taken very large, freezing out the dynamics of the D-branes wrapping it. From 
this point of view this was also studied recently in 

This paper is organized as follows. In Section |21 we discuss the matrix model approach 
and show in more detail the genus expansion for the different proposals. In Section IHl we 
study the 2-matrix model of type A2 and relate it to the flavor model. In Section 0] we work 
out explicitly the model with quadratic superpotential and compare with another approach 
using meson potentials. In Section El we conclude with same discussion. 



2. Flavors and Matrix Models 
2.1. Matrix Models With Flavors 

We consider an SU(A'c) gauge theory with Nf pairs of quark fields Qi and Qi, i = 1, ■ ■ ■ , A"/, 
in the fundamental and the anti-fundamental of SU(A'c) respectively. Note that the gauge 
superfield Wa contains a complex scalar $ in the adjoint representation of the color group. 
The Lagrangian density of this theory is given by 

L = Jd^9TT{Qie^Q, + Qie^Q,)+ J d^9Tr{W{<!>,Q,Q) + rWW^), (3) 
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with superpotential given by 



W{^, Q, Q) = W{^) + ^QQ - QmQ, (4) 

where the first term is a polynomial tree level superpotential for the adjoint Higgs field $ 
and m is the mass matrix for the fiavors. Classically the eigenvalues of the Higgs field and 
the quarks will be at critical points of this tree level superpotential. This breaks the gauge 
group to X J SU (Nj) where Nj is the number of eigenvalues of $ at the critical point I-th 
critical point. Quantum mechanically there will be a gluino condensate Sj for each of the 
factors. The goal is to find the effective superpotential for these condensates. 

The effective superpotential can be calculated according to the proposal of Dijkgraaf and 
Vafa, using matrix models. Here we take the superpotential of the gauge theory reduced to 
a point in supcrspace, replacing the number of colors by the matrix model size M of the 
number of ffavors Nf by Mf. So the fields are an M x M matrix $, an M x Mf matrix Q 
and a Mf x M matrix Q. The partition sum of the matrix model is given by 

z^e^ = /p$pgpge-^^^(*'«'^), (5) 

where W is the same function as above, but now of the matrix variables, and G is the 
unbroken gauge group. 

We denote by g the genus of the diagram, and by h the number of quark loops. By this 
we mean the number of boundary components coupling to the flavor matrices. Note that 
this is not the same as the number of holes (i.e. the number of holes coupling to the gauge 
multiplct). The dependence of the free energy of the parameters g^, M, and Mf can easily 
be extracted from the topology of the diagrams and can be written as an expansion in the 
genus and the number of quark loops as 

^ = JZas'^'-'M^J^.^^ig^M) = Y.9T\9sMffT,^n{gsM). (6) 

g,h g,h 

Let us therefore introduce parameters S = ggM and Sf = gsMf. 

We propose, analogous to DV, that the effective superpotential can be found from an 
expression of the form 

WefAS) = ^c^(^, ^/ = 0) + A^/^(5, ^/ = 0), (7) 

where is the planar contribution to the free energy for genus 0. This is a generalization 
of the pure gauge theory studied by Dijkgraaf and Vafa. 
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The planar contributions can be found from a large M expansion of the matrix model. 
All planar diagrams are summed by taking both the rank of the gauge group and the number 
of flavors to infinity, 

M, oo, ^ 0, S = gsM, Sj = gMf finite. (8) 

From the expansion above we see that this picks out the genus zero contribution. 

Let us relate the proposal above to the calculation of the matrix model in a large M limit 
taking Mf = Nf fixed, as in ^ and many papers following it. From the above expansion we 
find in the large M limit 

^ - g(^'=5,« - .iv,sr'^„,,.(5)). (9) 

Taking Sf = there are only two surviving terms: in the first sum we only get the J-'ofi con- 
tribution, while in the second sum we only retain J^o,i- Therefore the effective superpotential 
is given by 

f)7- f)7- 

W^ffiS) = N^^iS) + NfJ^oAS) = Nc^iS) + ^,=i(^). (10) 

We find the same expression as [6^ . The advantage of our proposal is that in keeping Mf = Nf 
finite, the contribution at x = 1 is subleading in Qs or equivalently 1/M. Therefore one 
cannot completely rely on large M expansion in the matrix model, which makes it harder to 
do actual computations, especially if one wants to consider general tree level superpotentials 
for the adjoint Higgs field. 

We should note that, as opposed to the proposal of [7, we do not relate Mf to Nf. 
Particularly, the ration Mf/M is not fixed in the large M limit. This is necessary, as we 
want to take 5*/ = at the end. Indeed, the matrix model quantities M and Mf are 
completely unrelated to the gauge theory quantities N and Nf. 



2.2. Loop Equations 

Let us first consider the matrix model with a large number of flavors, as described above. 
To calculate the matrix model integral, we integrate out the fundamentals. We are then left 
with an integral over $, 

Z = |p<De-^^^^(*\ (11) 

where 

y($) = H^($) + 5/log(<l>-m). (12) 
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We take the large M, Mf limit with gi^ — > such that the 't Hooft coupling S = QgM and 
Sf = QsMf is fixed. 

As is usual we introduce the resolvent defined as 

u;{x)^ ^Tr^—. (13) 
We can then derive the loop equation 

as a Schwinger-Dyson equation for the matrix model. Here we have to be careful that C is a 
cycle enclosing all the cuts but not the point x or m. The latter is because V'{z) has a pole 
at z — m. This loop equation can also be written in the form 

S{u;{xf) - VXx){u;{x)) + ^f{x) = 0, (15) 

where 

Note that if V were polynomial, this would be a polynomial in x. Also, in the large M limit 
we have (ci;(x)^) = {uj{x))'^ + 0{^). Therefore we can replace u!{x) by its expectation value. 
Introducing y{x) — —2Suj{x) + V\x) we find the spectral curve 

y^ = V\xf-f{x). (17) 

Unlike the usual matrix model however, the functions V'{x) and f{x) have a pole x — m. 
The pole term in f{x) is given by 

AS/^( I f Sf Sf \\\ ASSf , , 

— (Tr ^- > = ^-u{m). (18 

M\ \x-^\x-m ^-mJJ/ x-m 

As a result a multi-cut solution has the form 



y{x) = (^M{x) + ^JUi^^), (19) 

where the cuts lie along the intervals [x2i-i, X2i], M{x) is a polynomial with degree determined 
by the number of cuts and the degree of W, and 7 is determined by the pole at x = m to be 

7- ■ (20) 

yUiim - Xi) 
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Note that this coefficient is also such that the resolvent 

2Suj(x) = W'(x) + -^ y(x), (21) 

X — m 

has no pole at x = m. 

Rather than pursuing this direction, we want to proceed to a different method, relating 
the model to a special limit of a 2-matrix model. As it turns out several aspects of the 
solution can be seen much more clearly in this approach. 



3. Flavors From the 2-Matrix Model 

In this section we discuss the use of the 2-matrix model to calculate the effective superpoten- 
tial. We will interpret the two gauge group as the color and flavor symmetry group, which 
now are both gauged. To flnd the SU(A'"c) theory with fundamental flavors we take a limit 
forcing the flavor sector into a classical vacuum. 



3.1. The A2 Model 

We propose a description of the flavor model in terms of certain multi-matrix models. In 
particular, the A2 model. The idea is as follows. The A2 model is based on a quiver with 2 
nodes, both corresponding to a gauge group U(Aj). We have two bifundamentals, which we 
call Q = Q12 and Q = Q21. The tree level superpotential is given by 

Tn Wii^i) + Tt2 W2{^2) + Tri Q^iQ - Tra Q<^2Q- (22) 

Let us discuss the solution of this model. Part of the following discussion, and some more 
details, can also be found in j2]. We can write the matrix integrals in terms of eigenvalues 
Xij of the matrices We introduce the resolvents for the two adjoints, 

uj,{x) = ^Tii^— = (23) 

After integrating out the bifundamentals Q-,Qi we obtain [2] 

/n^V n (V-A,,,)^'^e-^S-.^^^(^-). (24) 

with Cij = 25ij — 6ij^i — 6ij^i the Cartan matrix of the A2 model. The forces, or equivalently 
equations of motion, can be written 

y,{x) = -gsdx,^,S{x) = W^ix) - a,SjU,{x), (25) 
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where Si = QgMi are the 't Hooft couphngs. It is convenient to introduce two scalar fields 
^Pi{x) on the x-plane related to the forces as d(pi{x) = yi{x)dx. The matrix integral can in 
fact be related to the rational CFT of type A2 in terms of these scalar fields. This CFT is 
known to have a structure of a W3 algebra formed by 2 currents W^'^^x) of spins s = 2,3. 
These currents also exist in the matrix model 

To write expressions for the currents, it is convenient to use a slightly different basis for 
the forces |2j, 

ai{x) = ti{x) — SiUi{x), 

a2{x) = t2{x) + SiUJi{x) - S2UJ2{x), (26) 

a-iix) = hix) + S2UJ2{x), 
where the ti are the classical contributions 

ti{x) = ^Wi{x) + ^Wiix), 

t,(x) = -^Wiix) + ^W^ix), (27) 

hix) = -Iwiix) - ^W^{x). 

These are not independent, but satisfy Y.iO'i = = Y^iti- The two spin s currents can then 
be given by jlHl 

n;{^)(a;) = ^zl)!^a,(x)^ s = 2,3. (28) 

^ 1=1 

They have a classical contribution yV^i\x) which can be found by replacing the with the 
ti. These currents satisfy loop equations of the form 

/^^(W(^)(.))=0. (29) 
Jc 2m X — z 

where C is a cycle enclosing all the eigenvalues but not the point x. The loop equations 
imply that 

{W^'\x)) - Wii\x) = w^'\x), (30) 

where w^'^\x) are regular functions of x. If the Wi{x) are polynomials of degree n + 1, one 
can show that w^'^^ is a polynomial of degree ns — n — 1. In the large M limit the expectation 
values factorize, so we do not have to distinguish between the resolvent as an operator and 
its expectation value in the expression for the current. 

The quantum spectral curve of this model can be written 

3 

Hiy - a,{x)) = 0. (31) 
1=1 
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Using the currents we can write this as a deformation of the classical curve — the above 
curve with the replaced by the tj — 

3 

Hiy - ti{x)) + yw^^\x) + w^^\x) = 0. (32) 

i=l 

This quantum curve can now be conveniently used to study the large M behavior of the 
model. 

r> sheet 1 




Figure 1 : The three sheets of the A2 curve, and the three cuts connecting them. The 
points P,Q,R are the points at regularized inhnity. 



3.2. Free Energy at Large M and the Effective Superpotential 

We will now study the large M limit by looking at the shape of the quantum curve. Let us 
first study the classical equations of motion of the system. In the following we will assume 
that there are no simultaneous solutions to Wl{\) = and H^2(-^) — 0- This is not a strong 
condition and will be sufficient for the rest of this paper. 

We can use the \J{Ni) x U(iV2) symmetry to diagonahze the matrices $j; their eigenvalues 
will be denoted Xij. For the corresponding blocks Qu and Qji in Q and Q the equations 
of motion give 

(Ai,7 - X2,j)Qij = 0, (Ai,, - X2,j)Qji = 0. (33) 
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This implies that the blocks Qu and Qji are nonzero only if Ai,/ = A2,j. We can therefore 
safely reduce the discussion to individual eigenvalues and henceforth suppress the extra index 
/. 

The remaining $j equations of motion are given by 

W[{\i) = -QQ, W^{X2) = QQ, (34) 

where the left-hand sides have to be interpreted as diagonal matrices. There are three types 
of solutions to these equations. 

When a block in QQ is zero, we find a solution for a critical point of the first node, 
W[{Xi) = 0. We will call these critical points eij. Similarly, we have solutions 62,7 corre- 
sponding to the critical points of the second node, pygl-^s) = 0. It remains to study the 
critical points when the eigenvalues are equal, Ai = A2 = A12. It follows from the equations 
of motion ()34j) and some elementary algebra that QQ = QQ} Therefore, these eigenvalues 
satisfy 

W[{\i2) + W!,{\r2) = 0. (35) 

The solutions to this equation will be called 612,7. 

The three different types of classical solutions can be related to the intersections of the 
three sheets comprising the classical curve. The intersection of the i-th and j-th sheet is 
given by the solutions to tj(x) — tj{x) = 0. 

Quantum mechanically, the eigenvalues for $j will spread around these classical solutions. 
In the large M limit they will form cuts in the complex x-plane. The contours around the 
cuts will be denoted Aij and ^12,7. The resulting quantum curve is sketched in Figure 
the 3 sheets correspond to the solutions y = aj(x). We also indicated the positions of the 
cuts around the different classical solutions. 

The number of eigenvalues, or rather the 't Hooft couplings Sij = ggMi j, can be measured 
by the contour integrals around the cuts, 

Si,i = (f y{x)dx, i = l,2; S^j = (f y{x)dx. (36) 

Integrating at infinity on the first and third sheet, we find that the total 't Hooft couplings 
Si = QsMi are given by the sums 

5*1 = ^ 5*1,/ + ^ S'12,7, S2 = '^82^ + '^Si2,i. (37) 
7 7 7 7 

^Assume, w.l.o.g., that W{{\i2) ^ 0. Then we can write, after a basis transformation, Q = [Qi 0) and 
Q — (^^\^ with Qi and Qi invertible matrices. We find QQ 7^ and, using (|34|) . it must be that Q — Qi. 
Therefore QQ = Q^^{QQ)Q = QQ, as QQ is central. 



9 



For our purpose we need to solve the 2-matrix model in the large M limit. At large 
M we can replace the sums over eigenvalues by the integral over eigenvalue distributions 
pi(Ai) and p2{^2) respectively. They can be found from the discontinuity in the correspond- 
ing resolvents. The large M limit of the free energy can be found from the saddle point 
approximation, 

-^o('S'i, = -Si J dXi pi{Xi)Wi{Xi) + ^CijSiSj J dXidXj pi{Xi)pj{Xj)log{Xi - Xj). (38) 

As in [1 , the free energy can be found from the curve by considering the change of the 
density when we remove an eigenvalue from the cuts. When we remove an eigenvalue from 
the cut at X = Cij, the density changes by 

S^AS^P^){k) = SS,j6{X, - 6,,,), t = 1,2. (39) 

Using this we find, for i = 1,2, a variation 

—— = -Wi{eij) + CijSj / dXj pj{Xj) log(Aj - 6^) = / {ai{x) - ai+i{x))dx. (40) 

As y{x) = ai{x) on the i-th sheet, we can write this as a contour integral over the noncompact 
curve Cij going through the appropriate cut, as indicated in Figure [H 

= / y{x)dx. (41) 

As the curves are noncompact we should regularize these integrals by ending the curves at 
finite points P, Q and R on the three sheets respectively. A similar discussion holds when 
we remove an eigenvalue at ei2,/. We now find the integral over ai{x) — a^{x), which can be 
written as an integral over the noncompact curve C12,/ connecting the first and third sheet, 

7^ = / yWx. (42) 

Let us now relate these to the effective superpotential of the quiver gauge theory. The 
original SU(A''i) x SU(iV2) gauge group will be broken by the distribution of the eigenvalues of 
the adjoint Higgs fields around the different classical vacua. The ranks therefore decompose 
analogously to 

^1 = E ^1,/ + E ^2 = E ^2,/ + E ^12,/, (43) 
II II 

analogous to (p?7jl . Here the SU(A''i2,/) factors are embedded diagonally in both original 
group factors. The effective potential for the A2 quiver theory should now be given in terms 
of the matrix theory quantities by 
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3.3. Decoupling the Flavor Group 

To relate the A2 model to the flavor model, we identify the group SU(A^2) with the flavor 
group. To decouple this group we need to take this sector in a classical vacuum. Doing so, 
the adjoint $2 of this node will play the role of the mass matrix m. To get the right masses 
we take the tree level superpotential of the flavor node such that 

^2(^2) =«^n(*2-m„). (45) 

a=l 

We will then take k, ^ 00. This ensures that the flavor node will not add any quantum 
corrections, so that indeed $2 will be constrained to the critical points rua. 

So we have to consider the limit k — > 00 of the cubic Riemann surface. In this limit the 
Riemann surface degenerates, and the third sheet will decouple. Note that for large k., all the 
cuts in sheet 3 corresponding to y — a3{x) are around the points W2 dominates. 

For large k around x — rUa, the resolvent S2U!2ix) will have the form 

S2UJ2{x) ~ ^Ka{x - m„) - ^^J kI{x - TUaY - 4:KaS2,a, (46) 

where is proportional to k. Ka and S2,a in the above equation will depend on x, but have 
finite values ai x = ma- Taking large k we have 

S2U;2{X) ~ l-Ka{x - m„) - l-^Kl{x-mar-^KaS2,a = + ^(-) " (^7) 

Therefore in the limit — > cxo we have 

S2u;2(x) = Y: (48) 

The above arguments would also allow extra regular terms, but the asymptotics u!2{x) ~ ^ 
for large x does not allow them. Therefore the above form is exact in the limit k, ^ oo. Of 
course this was expected; the eigenvalues of $2 are all in their classical vacua at A2 = m^. 

We see that in the hmit the branch cut dissapears, therefore the equation for the Riemann 
surface factorizes in this limit. We should then be interested in the part coming from the 
other two sheets. We have to change the coordinate y to absorb the large k part. The proper 
shifted coordinate on the Riemann surface can be given by 

y = 2y + a3{x). (49) 

We shift the ai_2 by the same amount, 

oi = 2ai + 03, 0,2 = 2a2 + 03. (50) 



11 



Note that we now have ai + a2 = 0. The equation [y — ai){y — 02) = for the decoupled 
first two sheets becomes 

{y - a,){y - a^) ^f- {W[{x) + S2u;2{x)f + f{x). (51) 

The quantum correction 

f{x) = -ASluj + 4SiUJiW[ + 4S1S2UJ1UJ2, (52) 

is a rational function with single poles at x — nia- This can be seen by writing this in terms 
of the quantum contribution to the spin 2 current, 

Slui{xf - S,S2u;i{x)u2{x) - S,uJi{x)Wi{x) = w^^\x) - S^u;2{xf + S2U2{x)W^{x) , (53) 

where w^'^^x) is a polynomial function of x. We then have 

where now C is a cycle containing all the eigenvalues of $1, but not the eigenvalues rUa and 
the point x. We then find that f{x) is regular away from x — rUa- As it only has at most a 
single power of a;2, it will have single poles &i x — ma- Near x — rUa it can be approximated 

f{x) ^ 4:SiS2,aUJl{ma) — ' • (55) 

X — nia 

In fact, this decoupled curve is precisely the one we found in the last section governing the 
matrix model with a large number of flavors. This shows that indeed as expected the two 
methods are equivalent. 

Let us change some notation appropriate for the present setting. We will denote Si = Sij, 
S'j'^ = 5*2,/ and Sj^ = S^j, and similar notations for the Nj. Let us now see what form the 
derivatives of the free energy reduce to. First we note that the CLl — (32 — ^(o-l — Ci2 

), which 

is half the difference of y on the two sheets. Therefore we can write the variation of the tree 
level superpotential with respect to the moduli 5*1,/ as 

dJ'o 1 



dSi 2 7c/^"^^"' ^'^^ 
where Cj is the noncompact curve in the decoupled double cover corresponding to Cij, 
running from Q to P. 

For the variation with respect to S^^, we note that the a2 — = |(a2 — 803). We have 
803 = W[ + 2W2 + 35'2u;2, but W2 is independent of S while the last term vanishes at 5*2 = 0. 
Therefore we can write 



/ (02(3;) — as{x))dx — - / {y{x) — W[{x))dx + constant. (57) 

JrUa 2 JrUa 
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We can safely ignore the extra constant for the calculation of the superpotential. In fact it 
is the contribution of the part of C2J lying on the decoupled third sheet, and represents the 
residue of the second node. Similarly, for the degenerate cuts connecting the first and the 
second sheet we find at 6*2 = 

{ai{x) — a^{x))dx = - f {y{x) — W[{x))dx + constant. (58) 



dSj^^ J ma 2 J ma 

The constants are related to the dynamics of the second node. As we are not interested in 
these, we can simply substract them. As they do not depend on the glueball expectation 
values this can be safely done. 

For the flavor model we are interested in the limit Sf = 0. The only relevant parameters 
are Sj = Sij. In the expressions above for the derivatives of the free energy we should 
substitute the solution for y at S2 = 0. This is simply the solution of the 1-matrix model 
with tree level superpotential Wi{x), and moduli determined by 5*/. 

The effective superpotential can now be calculated from with the proper interpreta- 
tion of the factors. Note that gives the decomposition of the ranks as 

= E + E Nla, Nf = Y: Nf^a- (59) 

I a a,± 

Hence we find that the different choices for = ±1 give rise to different unbroken flavor 
symmetries. 



3.4. Remarks on Seiberg Duality 

The A2 curve has a natural action of the A2 Weyl group 5*3. The part of this group compatible 
with our limit n 00 decoupling the third sheet is a Z2 symmetry interchanging sheets one 
and two. Looking at the action of this Z2 on the cuts they act on the 't Hooft couplings as 
Sj = —Sij, S2J = Si2j and S'12,7 = S2j- The curve after the Z2 action describes a similar 
gauge system but with different parameters. In terms of these dual variables the effective 
superpotential becomes 

Nj^ + Nt ^ + N7 ^ = -N,^ + N7 ^ + Ni ^ (60) 



We read off from this that this can be interpreted as the effective superpotential for a dual 
theory with Nj = —Nj while N^^ = Nj^. This implies 

Nc = J:Ni + E^M = -E^/ + E^i = ^/ - Nc, Nf = Nj. (61) 
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Note that for this to have a gauge theory interpretation we need Nj > Nc, as the number 
of colors is always positive. This is exactly the transformation of Seiberg's electric-magnetic 
duality |i4Sj. Indeed this Weyl group symmetry has been connected to Seiberg-like duality 
in the context of D-branes on Calabi-Yau geometries with ADE geometry in ^\ . 

Note that the quantum curve has two types of degenerate cuts connecting the third "fla- 
vor" sheet, corresponding to = ±1. These have a different physical interpretation. This 
can be seen form the behavior of the Nc dependence. The duality described above inter- 
changes these two cuts. In the Seiberg duality, the magnetic model also is slightly different, 
as it has extra gauge invariant fields, and an extra term in the tree level superpotential. 
Furthermore the original Seiberg duality has massless quarks, so we need to take a limit 
m ^ 0. 



4. The Gaussian Model 
4.1. Matrix Model Calculation 

Let us for specifity consider the simplest example of U{Nc) super Yang-Mills with Nj massive 
flavors and an adjoint, with a gaussian tree level superpotential for the adjoint. That is, we 
take 

iyi($i) = 1$?. (62) 

The superpotential W2 of the second node will be taken in terms of the masses for the 
flavors as above. 

The tree level superpotential has only a single vacuum at $1 = 0. Hence there will be 
only a single glueball Sij = S. For 6*2 = the solution for y{x) is has the well known half 
circle form, 

y{x) = ^Wi{xy-f{x) = ^g^x^ - 4gS. (63) 

(We leave out the tilde as we are only considering the decoupled system). 

We take the noncompact cycle C = Cj, connecting the point running from Q to P passing 
through the 5* cut. We will take a regularization where P and Q are taken at A. We then 
have, at S2 = 0, 5*1 = S, 

§4/^.(xMx^-Slog-fj + 5. (64) 
For the derivative with respect to Sf^ at S*/ = we have to take the integral from the point 



14 



nia to either P or Q. We find, up to a S* independent constant, 



9St 2 



{y{x) -W[{x))dx (65) 



ma\ 1 



S'logf^) - ^5 + S'log7+(aa5') + 



A J 2 -'-^v » / ' 27+(«a^)^ 



where a„ = and we have introduced the functions 

1 1 



7±(x) = - ± -VT^. (66) 
Similarly, we can calculate the derivative with respect to Sj^^, for which we find 

-f^ = n {y{x)-W[{x))dx (67) 

UO j g, Z J ma 



\gmaAJ 2 2-f_{aaS)' 

From our proposal ((Tj), we can give the following expression for the effective superpotential, 
noting that Ni = Nc - Nj due to 

W^jfiS) = -iV,^(log^ - l) + E ^(log7.„(aa^) + ^^^^g^ - 1). (68) 

where ea = ±1 according to whether the integral was taken to Q or P, and the scale A is 
given in terms of the cutoff as A^^'^ = K^^^'^f g'^" det m. This is precisely the form found by 
an explicit summation of the planar diagrams with up to one boundary |Hl I15j. 

The actual effective superpotential can be found by minimizing the above with respect 
to 5*1, dsWeff = 0. The critical points are determined by the solutions to the equation 



g \ N, ^ 



A3 



UleSc^aS). (69) 



a=l 



At a critical point the logarithmic terms in the effective superpotential cancel, and it can be 
written 

i^w-^^s+i:^^^. (TO) 

where S" is a solution to the equation above. 
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For simplicity we will now consider the case where all masses are equal to m and all 
= 1. This gives an equation for S which can be written 

2Nc-Nf Nc-Nf 

Here we used that 5* = can not be a solution. The expression for the effective superpotential 
at the critical point is given by 

A3 / C / a3x -V^\ 

Weff,crU=Y[(^N^-Nf)j^ + Nf[Y) ' j' (72) 

where for S we should insert a solution to the above equation. 



4.2. Calculation Via Meson Potential 

Another way to calculate the answer is to integrate out $ and express everything in terms of 
the meson field X = QQ. This procedure was used in |6] for the case of a single flavor, and 
in for more flavors. It is useful though for us to review it as it will give an interpretation 
of the two flavor cuts in terms of meson expectation values. 

After integrating out $, we obtain a tree level superpotential which is quadratic in terms 
of the meson field X. To get the full effective superpotential, we need to include the non- 
perturbative Affleck-Dine-Seiberg potential jlU]. We then get an effective superpotential for 
the meson field of the form 

fA^^^-^f\w^ 1 



W,„^iK-N,)[^^^) ■ +Tr,nX--TrX\ (73) 

The equation for a critical point of this superpotential is 



l-mX + -X' = 0. (74) 
det X J g 

At the critical point we can write the superpotential in a slightly simpler form, 

2N^- Nf /A^^-~^f\N^ 1 , , 

= '+2T.,„X (75) 

More generally, we can diagonalize the meson matrix X. We can actually assume that 
X is diagonalized simultaneously with the mass matrix m. This follows from the equations 
of motion for the flavors, 

^ = ^Q-Q^ = ^^-^^- (76) 
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Note that these equations of motions are exact quantum mechanically, as the flavors appear 
only quadratic in the action. We will therefore call the eigenvalues of the meson field X^. 
The equations for the critical points for the eigenvalues can now be written 

a-maXa + -Xl = Q, (77) 



where we have introduced a = , y . We can now solve for the eigenvalues of the 



t3Nc-Nf \ Nr-N 
detX 

meson field in terms of a, finding 

Xa = T, (78) 

"^a7Ea(«aCr) 

where = ±. Due to the obvious relation detX = Ha^a, o" has to satisfy the relation 

n 7e.(«a^) = ;^3^ ^ ■ ■ (79) 

We now find that if we relate the scales as A^^'^"^-'' = and we replace a by S, this 

is precisely the equation ()69|) for 5*. Moreover with this identification also the effective 
superpotentials (fTHjl and (f7H|) agree. 

As the choice of is related to the choice of sheet on which the point corresponding to 
rUa, we see that the two different branch cuts in the A2 model near x = rua correspond to 
the choice of eigenvalues for the meson field. 



5. Conclusions and Discussion 

We studied effective glueball superpotentials for A/" = 1 super Yang-Mills theories with 
fundamental matter and an adjoint Higgs field. We calculated these using matrix models 
at large N with a large number of flavors. In particular, we used the 2-matrix model in a 
special limit. As noted earlier, this allows us to deal exactly, and purely in a large N limit, 
with arbitrary tree level superpotentials for the adjoint field. Although we only worked out 
in detail the Gaussian case, more general superpotentials are straightforwardly added, but 
technically more difficult to solve exactly. 

The description of flavors using the A2 quiver matrix model has some intriguing proper- 
ties. The presence of two types of cuts in the A2 model connecting the sheet corresponding 
to the flavor As noted above, we can see a sign of Seiberg-like duality from the Z2 action 
interchanging two sheets. It would be interesting to pursue this further. 
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